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SELF— S HIELDED GRADIENT FIELD COIL FOR MAGNETIC RESONANCE IMAGING 

The following relates to the diagnostic imaging arts. 
It particularly relates to magnetic field gradient coils for 
use in short -bore magnetic resonance imaging scanners, and will 
be described with particular reference thereto. However, the 
5 following relates more generally to magnetic field gradient 
coils for various types of magnetic resonance imaging scanners, 
and to dedicated gradient coils for imaging reduced volumes of 
interest . 

In a magnetic resonance imaging system, spatial 

10 encoding of the magnetic resonances is typically accomplished 
by generating magnetic field gradients in an imaging region. 
These gradients are produced by a set of magnetic field 
gradient coils. Gradient coils are typically provided for 
producing independent magnetic field gradients in each of the 

15 transverse x- and y-directions , and in the axial z-direction. 

For various applications, such as producing gradients 
in short -bore magnetic resonance imaging scanners, or producing 
dedicated localized gradients over a reduced field of view 
(FoV) for imaging specific organs, it is desirable to employ 

20 relatively shorter gradient coils. In scaling gradient coils to 
reduced length while retaining high gradient peak amplitudes 
and linearity and bore diameter, coil turns become increasingly 
compressed near the coil extremities. Moreover, imaging 
distortions and artifacts due to eddy currents induced in a 

25 surrounding cold shield become more problematic. Furthermore, 
adaptive or variable field of view gradient coils are desired 
to adapt a given magnetic resonance scanner or magnetic field 
gradient coil set for a wider range of applications. 

The present invention contemplates an improved 

30 apparatus and method that overcomes the aforementioned 
limitations and others. 
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According to one aspect, a gradient coil is disclosed 
for a magnetic resonance imaging apparatus. A primary coil 
defines an inner cylindrical surface. A shield coil defines an 
outer cylindrical surface coaxial ly aligned with the inner 
5 cylindrical surface and having a larger cylindrical radius than 
the inner cylindrical surface. A plurality of coil jumps 
electrically connect the primary and shield coils. The coil 
jumps define a non-planar current -sharing surface extending 
between an inner contour coinciding with the inner cylindrical 

10 surface and an outer contour coinciding with the outer 
cylindrical surface. The primary coil, shield coil , and coil 
jumps cooperatively define a current path that passes across 
the current -sharing surface between the inner and outer 
contours a plurality of times. 

is According to another aspect, a method is provided for 

producing a magnetic field gradient in a magnet bore of a 
magnetic resonance imaging apparatus. An electrical current is 
circulated through a primary coil that defines an inner 
cylindrical surface. The electrical current is circulated 

20 through a shield coil that defines an outer cylindrical surface 
coaxially aligned with the inner cylindrical surface and having 
a larger cylindrical radius than the inner cylindrical surface. 
The electrical current is communicated between the primary and 
shield coils via a plurality of non-planar coil jumps a 

25 plurality of times. 

According to yet another aspect, a method is provided 
for designing a gradient coil for a magnetic resonance imaging 
apparatus. Current densities are computed on coaxial inner and 
outer cylindrical surfaces using constraints on the stored 

30 energy and gradient field linearity and uniformity. The 
computed current densities are generally non-zero at 
non-coplanar boundary contours that define coil edges. Primary 
coil turns, shield coil turns, and coil jumps are selected to 
approximate the computed current densities on the inner and 

35 outer cylindrical surfaces. 
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One advantage of designing the gradient coil with 
shared turns on the primary and shield coil is that it allows 
shorter shielded gradient coil designs. 

Another advantage resides in reduced magnetic field 
5 gradient distortion due to eddy currents in the cold shield or 
other conductive elements of the magnetic resonance imaging 
scanner. 

Another advantage resides in reduced current winding 
densities at the extremities of the gradient coil . 
10 Yet another advantage resides in reduced I 2 R heating, 

particularly in regions of high winding density. 

Numerous additional advantages and benefits will 
become apparent to those of ordinary skill in the art upon 
reading the following detailed description of the preferred 
15 embodiments. 



The invention may take form in various components and 
arrangements of components, and in various process operations 
20 and arrangements of process operations. The drawings are only 
for the purpose of illustrating preferred embodiments and are 
not to be construed as limiting the invention. 

FIGURE 1 shows an exemplary magnetic resonance 
imaging apparatus, in partial cross-section, and auxiliary 
25 components for performing magnetic resonance imaging. 

FIGURE 2 shows an enlarged cross-sectional view of 
the magnetic resonance imaging apparatus which shows a flared 
arrangement of primary and shield gradient coils and a 
current-sharing surface for passing current back and forth 
30 between the primary and shield gradient coils. 

FIGURE 3 diagrammatically shows a sectional view of 
the geometrical arrangement of the primary and shield gradient 
coils of FIGURE 2, with a suitable dimensional nomenclature 
superimposed. 
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FIGURE 4 diagrammatically shows a perspective view of 
the geometrical arrangement of the primary and shield gradient 
coils . 

FIGURE 5 shows an unfolded half-quadrant layout of 
5 the X-gradient primary and shield gradient coils and the 
interconnecting coil jumps. 

FIGURE 6 shows exemplary isolated primary and shield 
coil turns of the X-gradient coil of FIGURE 5. 

FIGURE 7 shows exemplary communicating primary and 
10 shield coil turns of the X-gradient coil of FIGURE 5. 

FIGURE 8A shows a plot of gradient coil system energy 
versus the z-component of the current density at the edge of 
the primary coil for an exemplary transverse gradient coil 
design which has G x = 30mT/m and the eddy current effect is - 
15 0.22%. 

FIGURE 8B shows a plot of gradient coil system energy 
versus residual eddy current effect for an exemplary transverse 
gradient coil design which has G x =30 mT/m and J 2p , e dge=17 kA/m. 

FIGURE 9 shows a plot of the z-component of the 
20 current density for the primary and shield coils layout shown 
in FIGURE 5. 

FIGURE 10 shows a plot of X-gradient non-uniformity 
at a radius of 0.001 meters for the primary and shield coils 
layout shown in FIGURE 5. 
25 FIGURE 11 shows a plot of X-gradient non-linearity at 

z=0 meters as a function of radial position for the primary and 
shield coils layout shown in FIGURE 5. 

FIGURE 12 shows a plot of residual eddy current 
distribution as a function of axial position Z produced by the 
30 coil layout shown in FIGURE 5. 

FIGURE 13 shows a plot comparing eddy current effect 
as a function of axial position Z for the continuous current of 
FIGURE 9 and the discretized coils of FIGURE 5. 
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FIGURE 14 shows an unfolded half-quadrant layout of a 
coil P2 formed by selecting fourteen turns on the primary coil 
of FIGURE 5 that do have jumps. 

FIGURE 15 shows an unfolded half-quadrant layout of a 
5 coil Pi formed by disconnecting fourteen isolated primary coil 
turns of the primary gradient coil P2 of FIGURE 14. 

FIGURE 16 shows an unfolded half -quadrant layout of a 
second shield coil S2 that operates in conjunction with the 
coil PI to define a second field of view. 
io FIGURE 17 shows a plot of X-gradient non-uniformity 

for the two discrete fields of view of the gradient coil 
illustrated in FIGURES 5 and 14-16. 

FIGURE 18 shows a plot of X-gradient non-linearity 
for the two discrete fields of view of the gradient coil 
15 illustrated in FIGURES 5 and 14-16. 

FIGURE 19 shows an unfolded half -quadrant layout of a 
second shield coil S2 ' that operates in conjunction with the 
coils of FIGURES 5, 14, and 15 to define a discretely variable 
gradient coil with three fields of view. 
20 FIGURE 20 shows two gradient coils driven in 

parallel, with nomenclature specified for computing a slew rate 
for the parallel-combined coils. 

FIGURE 21 shows a suitable configuration of two 
primary coils located on a single surface and two shield coils 
25 designed as described herein that provide a continuously 
variable field of view. 

FIGURE 22 shows a plot of the z-component of the 
current density for primary and shield coils providing a 32cm 
DSV field of view. 
30 FIGURE 23 shows a plot of the z-component of the 

current density for primary and shield coils of a correction or 
shim coil that cooperates with the 32cm DSV coil of FIGURE 22 
to produce a continuously variable field of view. 
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FIGURE 24 shows an unfolded half-quadrant layout of 
the primary and shield gradient coils and the interconnecting 
coil jumps of the discretized current densities of the 32cm DSV 
coil of FIGURE 22. 
5 FIGURE 25 shows an unfolded half-quadrant layout of 

the primary gradient coil of the correction coil. 

FIGURE 26 shows an unfolded half -quadrant layout of 
the shield gradient coil of the correction coil. 

10 

With reference to FIGURE 1, a magnetic resonance 
imaging scanner 10 includes a cylindrical main magnet 12, which 
is preferably superconducting and cryoshielded. The main magnet 
12 defines a magnet bore 14 inside of which a patient or other 

15 imaging subject is placed for imaging. The main magnet 12 

produces a spatially and temporally constant and uniform main 
magnetic field oriented along a longitudinal axis of the bore 
14. Instead of a superconducting magnet, a non- superconducting 

magnet can be used. Moreover, a vertical magnet, an open 
20 magnet, or other type of main magnet can be employed instead of 
the illustrated horizontal cylindrical main magnet 12. In one 

preferred embodiment, the main magnet 12 is a short-bore 

magnet . 

Magnetic field gradient coils include a primary 
25 gradient coil 16 and a shield gradient coil 18 that 

cooperatively produce magnetic field gradients in the bore 14 
for spatially encoding magnetic resonance signals, for 
producing magnetization-spoiling field gradients, or the like. 
Preferably, the magnetic field gradient coils 16/ 18 include 
30 coils configured to produce magnetic field gradients in three 
orthogonal directions including transverse x- and y-directions . 
In addition to the shield coil 18, a cold shield 20 provides a 

high conductivity eddy current surface for residual gradient 
fields thus protecting the magnet coils still further away. In 
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a suitable embodiment, the cold shield 20 is integrated into a 
housing of the main magnet 12. 

A whole body radio frequency coil assembly 22 
generates radio frequency pulses for exciting magnetic 
resonances. The radio frequency coil assembly 22 also serves to 
detect magnetic resonance signals. Optionally, additional local 
radio frequency coils or phased radio frequency coil arrays 
(not shown) are included for exciting and/or detecting magnetic 
resonances at localized areas in the bore 14. 

Gradient pulse amplifiers 30 deliver controlled 
electrical currents to the magnetic field gradient coils 16 , 18 
to produce selected magnetic field gradients. A radio frequency 
transmitter 34, preferably digital,, applies radio frequency 
pulses or pulse packets to the radio frequency coil assembly 22 
to generate selected magnetic resonance excitations. A radio 
frequency receiver 36 also coupled to the radio frequency coil 
assembly 22 receives magnetic resonance signals. If more than 
one radio frequency coil is provided (such as a local coil or 
phased coil array) , then different coils are optionally used 
for the magnetic resonance excitation and detection operations. 

To acquire magnetic resonance imaging data of a 
subject, the subject is placed inside the magnet bore 14, 
preferably at or near an isocenter of the main magnetic field. 
A sequence controller 40 communicates with the gradient 
amplifiers 30 and the radio frequency transmitter 34 to produce 
selected transient or steady state magnetic resonance 
configurations in the subject, to spatially encode such 
magnetic resonances, to selectively spoil magnetic resonances, 
or otherwise generate selected magnetic resonance signals 
characteristic of the subject. The generated magnetic resonance 
signals are detected by the radio frequency receiver 36, and 
stored in a k- space memory 42. The imaging data is 
reconstructed by a reconstruction processor 44 to produce an 
image representation that is stored in an image memory 46. In 



WO 2004/046745 



PCT/IB2003/004802 



one suitable embodiment the reconstruction processor 44 
performs an inverse Fourier transform reconstruction. 

The resultant image representation is processed by a 
video processor 48 and displayed on a user interface 50, which 

5 is preferably a personal computer, workstation, or other type 
of computer. Rather than producing a video image, the image 
representation can be processed by a printer driver and 
printed, transmitted over a computer network or the Internet, 
or the like. Preferably, the user interface 50 also allows a 
10 radiologist or other operator to communicate with the magnetic 
resonance sequence controller 40 to select magnetic resonance 

imaging sequences, modify imaging sequences, execute imaging 
sequences, and so forth. 

With reference to FIGURE 2-4, the magnet 12 of the 
is magnetic resonance imaging scanner 10 is preferably a 
short-bore magnet defining a cylindrical magnet bore 14. The 
primary and shield gradient coils 16, 18 coincide with and 
define axially concentric cylindrical surfaces 60, 62, 

respectively. The inner and outer concentric cylindrical 
20 surfaces 60, 62 have radii Rp, R s , respectively. The shield 

gradient coil 18 extends axially a half-length L s along the 
outer cylindrical surface 62, while the primary gradient coil 
16 extends axially a half-length 1^ along the inner cylindrical 
surface 60. The cold shield 20 has a length 2L CS . Preferably, 
25 the gradient coils 16 , 18 and the cold shield 20 are 
bilaterally symmetric along the axial direction, and L s , L p , Ii cs 
designate half-lengths. 

To provide reduced length, increased openness, and 
improved shielding, the axial half-length L s of the shield 

30 gradient coil 18 is preferably greater than the axial half- 
length L p of the primary gradient coil 16. The two gradient 
coils 16, 18 are connected along a non-planar current -sharing 
surface 64 that extends between an inner contour 66 coinciding 
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with the inner cylindrical surface 60 , and an outer contour 68 
coinciding with the outer cylindrical surface 62. In a 
preferred embodiment, the non-planar current -sharing surface 
corresponds to a curved surface of a frustum of a cone, where 
5 the cone has an angle a respective to an axis of the axially 
concentric cylindrical surfaces 60, 62. The angle a corresponds 

to an arctangent of a ratio of a difference between the 
cylindrical surface radii R^, R s and a difference between the 

primary and secondary or shield gradient coil half-lengths L p , 

io L s . However, the non-planar current -sharing surface can also 

have a parabolic, piece-wise conic, or other shape. 

With continuing reference to FIGURES 2-4 and with 
further reference to FIGURE 5, the primary gradient coil 16 

includes a fingerprint -like pattern of conductive windings or 
15 current paths coinciding with the inner cylindrical surface 60, 

while the shield gradient coil 18 includes a fingerprint -like 

pattern of conductive windings or current paths coinciding with 
the outer cylindrical surface 62. 

With continuing reference to FIGURES 2-5 and with 
20 further reference to FIGURES 6 and 7, the primary gradient coil 
16 includes isolated primary coil turns 70 that connect with 
other primary gradient coil turns on the surface 60. The 
primary gradient coil 16 also includes communicating primary 
coil turns 72 that each connect with at least one coil jump 74 
25 that coincides with the current - sharing surface 64. Similarly, 
the shield gradient coil 18 includes isolated shield coil turns 
80 that connect with other shield coil turns on the surface 62, 
as well as communicating shield coil turns 84 that each connect 
with at least one coil jump 74. 

30 As shown in FIGURE 6, each isolated primary or shield 

coil turn 70, 80 defines a current loop on the cylindrical 
surface 60, 62, respectively. As shown in FIGURE 7, the 
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communi eating primary and shield coil turns 72 , 82 and two coil 
jumps 74 define a current loop that passes back and forth 
across the current -sharing surface 64 via the two coil jumps 
74. 

5 Those skilled in the art will recognize that the 

gradient coils 16, 18 have substantial advantages over past 
coils, especially for short-magnet designs. By connecting the 
primary coil 16 and secondary or shield coil 18 with a 
plurality of coil jumps 74, high coil turn densities near the 
10 coil edges are reduced. Moreover, by axially flaring the shield 
coil 18 relative to the primary coil 16 improved shielding is 
obtained versus a design using primary and shield coils of 
equal lengths. The axial flaring of the shield coil 18 is 

accommodated by the non-planar current -sharing surface 64 to 

15 permit a combination of a flared shield and multiple coil 
jumps . 

With returning reference to FIGURE 3, a suitable 
method for designing the gradient coils 16, 18 uses the 

following parameters. The angle a is called a flare angle. The 
20 axial coil half-lengths L p/ L s , and the cold shield half-length 
L cs are shown and were discussed previously. Preferably, 
L p <L s <L cs , with L cg substantially larger than L s . The primary coil 
16 coincides with and defines the inner cylindrical surface 60 
which has a radius Rp. The shield coil 18 coincides with and 
25 defines the outer cylindrical surface 62 which has a radius R s . 
The cold shield 20 is preferably cylindrical with a radius R^, 
and is preferably coaxial with the cylindrical surfaces 60, 62. 
For a preferred embodiment, Rp<R B <R cs . 

The gradient coil design goals include producing: (i) 
30 given gradient strength G x or G Y at the isocenter of the 
magnet bore 14; (ii) high slew rate; (iii) good quality 
gradient magnetic field characteristics, such as linearity and 
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uniformity, within a given diameter spherical volume (DSV) ; and 
(iv) good shielding. It is recognized that leakage of the 
magnetic field into the conductive cold shield 20 will be 

present. This leakage is responsible for the residual eddy 
5 currents induced in the cold shield 20 when the gradient coil 

16, 18 is operated. These residual eddy currents can affect the 
quality of the gradient field inside the DSV and are 
characterized by the residual eddy current effect, which is 
defined as follows: 

o 

-peddy current/ \ 

ed_ef(r)| DSV = ^ Q " ra , ; rj ] DSV *100% (1) 



where the symbol | DSV indicates that the residual eddy current 

effect is evaluated at the surface of the DSV. One skilled in 
15 the art knows that if the above characteristic is reasonably 
constant as a function of spatial position, r, then the 

residual eddy current effect can be substantially corrected 
with gradient pre-emphasis . Thus, the absolute level of eddy 
current magnitude, although generally desired or preferred to 
20 be as small as possible, is not necessarily the most important 
aspect . 

In designing a self -shielded gradient coil, which 
from the beginning has finite lengths of the primary and shield 
coils, shielding conditions for the normal component of the 

25 total magnetic field B r (r) on the surface of the shield coil are 
typically imposed according to: 

*,(r)U=0 

s(k) = F s (k) + R(k)F p {k)-+0, R(k) = Y\^f p) (2) 

30 F P Ak)= \f^'\z)cos(kz)dz 

-V 
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where f£ PtS) (z) is the ^-component of the current density on the 
primary (p) or shield (s) coil, L ps are the radii and 

half-length of the primary (p) or shield (s) coil, 
respectively, and l[(kR) is a derivative of a Bessel function of 
5 the first kind i\ (fcff) . Equation (2) holds for any k, and is 
substantially restrictive. The shielding error function e(k) is 
a measure of how well the gradient coil is shielded. To design 
a well -shielded gradient coil Equation (2) leads to high 
current density on the primary coil with a negative lobe near 

10 its edge, leading therefore to high inductance. 

If the residual eddy current effect is non-zero it 
should have small variation over the imaging volume to be 
correctable. This means, that the variation of the residual 
eddy current effect should be small, even if the residual eddy 

15 current effect is about, say, 10%. Therefore, the coil is 
designed to satisfy the conditions that there is no high 
current density near the ends of the coils and the residual 
eddy current effect is correctable. 

The gradient coil design includes the following 

20 process operations: (i) allowing the current to flow off the 
end of the primary coil 16 to the shield coil 18 along the 
conical or other non-planar current - sharing surface 64 disposed 
therebetween; (ii) the cold shield 20 is included in the coil 
design; and (iii) the residual eddy current effect is 

25 considered as a variable in the design together with gradient 
field quality characteristics. The latter will be introduced as 
point constraints. 

In the following, an exemplary X-gradient coil design 
is described. Those skilled in the art can readily adapt the 

30 exemplary X-gradient coil design to the designing of a 
Y-gradient coil. The current distribution on the X-gradient 
coil in terms of cylindrical coordinates has the following 
form: 
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J(r) = e p J p (r)+e^(r) + e s y 2 (r) 

J p (r) = J™fr){S(z-Z(j?))+S(z + Z<j?))}®(p-R p )®(R s -/>)sin(p) 
J„(r) = {f™(z)S(p-R p )®(z + L p )®(L p -z)+f^(z)S(p-R,)®Cz+L 1 )®iL 1 -z)}cos<» 
•/,(*) = {f< p \z)S(p-R p )®iz+L p )®(L p -z) + fy( Z )5(p-R,)®(z + L s )®(L, -z)jsin( <P ) + 
J™(p){5(z-Z(j>y)-S(z + Z(p))}®(p-R p )®(R, -yo)sin(9>) 

Z(p) = L p +(p-R p )com, co\{&) = ^^- 

k ,~ k p 

(3) 

The continuity equation VJ(r) = 0 applied to the current J(r) 
given in Equation (3) above leads to the following equations: 



ggg^ .o- >J o. M .£, 

^C*,) = / 2 (p) (i P ), J 00 ^)- (4) 

C=R p f™(L p ) = -R t f?HL p ) 



where the functions //**(z) are suitably represented as Fourier 
series according to: 

/,<»(*) = 2>« sin(^z), /«(z) - sin(^-z) (5 ) 

with the coefficients a^ s) determined through extremization 

of the following functional : 

W = E -^(Bf (r, ) - B,)A, - £(2ff (r,) - C,^ ( r ,))A y - 

/ i 

J? (6) 
(fJ>(L p )-A p )A p -(f;( L ,)+^-A p )A, 
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where: E represents total energy stored in the magnetic fields 
of the gradient coil 16, 18 and the cold shield 20; 5 z C0l7 (r.) is 
the z -component of the coil magnetic field at a set of points 
r,. inside the DSV where the field has predetermined values B t ; 

the z-component of the magnetic field produced by the 
cold shield 20 at a set of points r y on the surface of the DSV; 

is the z-component of the coil's magnetic field at the 
same set of points; and C j is a set of numbers that govern the 
residual eddy current effect. Finally, X t and A y are the 
corresponding Lagrange multipliers. The quantity A is a 
desired value of the z-component of the current density on the 

... K 

primary coil at its edge, and the quantity Z-A p is the value 

of the z-component of the current density at the edge of the 
shield coil. The continuity conditions are represented by the 
last two terms in Equation (6) and A p , A s are the 
corresponding Lagrange multipliers for the current continuity. 

The quantities A p and C y are not variables and can be chosen 
arbitrarily. In deriving the Equation (6) , we have assumed 
that the normal component of the total magnetic field vanishes 
on the surface of the magnet cold shield 20, which has radius 

R cs , which is assumed much longer than the length of the 
gradient set. This gives the following equation for the cold 
shield residual eddy current: 
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f™(z) = -)F™(k)cos(kz)dk 
71 o 

F^ (&) = —R{k, R p , (Jc) — R(k, R s » R^ * (Jc) 

^MlM^ri^ dp (7) 
Res i/mcs) p 

Rcsli (Was) 

where the Fourier transforms F ps (k) of the primary and the 
shield current distributions are given in Equation (2) . 
5 Variation of the functional W of Equation (6) with respect to 
a%, a 3 n determines the solution in terms of continuous current 
distributions . 

Optimization of the functional given in Equation (6) 

with respect to a { n p) and a ( n s) determines a suitable solution by 
10 adjusting the value of the current density at the edge of the 
primary coil . After a solution is found the continuous current 
distributions on primary coil, shield coils and the 
interconnecting surface are discretized such that each loop on 
the primary coil carries the same current I, and each loop on 
15 the shield coil carries opposite current. The discretized coil 
is tested for its performance: gradient strength at isocenter, 
linearity and uniformity of the gradient strength over the DSV, 
and slew rate and residual eddy current effect. 

The continuous current distributions are discretized 
20 and approximated by coils 70, 72, 74, 80, 82 carrying current 
along a current path. A suitable fingerprint pattern of the 
primary and shield coils 16, 18 and coil jumps 74 determined by 
the forgoing process is shown in FIGURE 5. The left and right 
half of the coil are the primary and shield coils 16, 18 

25 respectively joined by coil jumps 74 at fixed angular positions 
(J) along the current -sharing surface 64 that connects the 
concentric inner and outer cylindrical surfaces 60, 62. 
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With continuing reference to FIGURE 5 and with 
further reference to FIGURES 6 and 7, the two kinds of current 
paths, namely the isolated primary and shield coil turns 70, 80 
exemplarily shown in FIGURE 6 and the communicating primary and 
5 shield coil turn 72, 82 shown in FIGURE 7 are generated. The 
isolated current loops 70, 80 reside only on the primary or 
shield coil 16, 18, that is, only on one of the two cylindrical 
surfaces 60, 62. The communicating current paths 72, 82 define 
current loops that reside on both the primary and shield coils 

10 16, 18, that is, on both of the two cylindrical surfaces 60, 
62, with coil jumps 74 defining current jumps or 
interconnections between them. 

The B z field from all of these paths is calculated to 
determine residual eddy currents in the cold shield 20, which 

15 can be compared to the continuous current residual eddy current 
effects on the DSV. Each isolated loop 70, 80 resides on 
either the primary coil 16 or shield coil 18 and is closed on 
itself. Examples of these loops are shown in FIGURE 6, which 
shows an isolated primary loop 70 on the inner cylindrical 

20 surface 60 and an isolated shield coil loop 80 located in the 
first quadrant, for example. Each loop carries a current I. The 
current on the isolated primary coil- loop 70 flows clockwise, 
while the current on the shield coil loop 80 flows 
counterclockwise. Exemplary communicating primary and shield 

25 loops 72, 82 defining a current loop passing back and forth 

across the non-planar current sharing surface 64 is shown in 
FIGURE 7. 

A suitable process for calculating the magnetic field 
and associated residual eddy current in the magnet cold shield 
30 20 from the discretized coil 16, 18 is the path- integral 
approach. The coil is discretized, so an equation describing 
each loop is known, i.e. the equations q>-f p {z\ <P = f s (z). For 
the isolated loops shown in FIGUER 6, the conditions 
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f P CP 1 ) = f P Cp2) = 0, f s (si) = f s (s2) = 0 apply. For the communicating 
loops shown in FIGURE 7, different conditions 
/,Cpl) = /,C*l) = 0, f p (L p ) = <p f , f s (L s ) = <p f apply, where <p f is the 
angle at which the current -sharing coil jump occurs. Rather 
than the parameters p x , p 2 , s x , s 2 , the description henceforth 
uses z i and z f . 

The current distribution produced by the 
communicating loops 72, 82 shown in FIGURE 7 is given by: 

J(r) = e p J p (r) + e, J, (r) + e z J z (r) 

J*(r) = Q p R p f p (z)^ 

+ 0 S R s f s (z){8{<p - f s (z)) + 8{<p + f s (z))}8(p - i? 5 )©(z> - z)©(z - z/ ), ( 8 ) 
•A (r) = {<%> - f p (z)) - <%> + f p (z))}d(p - i? p )©(z; - z)©(z - z* ) 

+ B 5 {5{<p - (z)) - <?(p + /, (z))}S(p - ^ )©(z} - z)©(z - z? ) 

+^(/7)^(z-Z(p))0( / 7-i? p )©(^ - p ){s(<p-<p f )-.S(<p + <p f )} 

where in Equation (8) Q p ,B p ,Q s ,B s are unknown constants and 

J p (p),J z (p) are unknown functions to be determined from the 
continuity equation VJ(r) = 0. Thus, 

J z (p) = J p (p)cot(S), J p (p) = - (9) . 

P 

This allows expressions for all components of the vector 
potentials, and therefore for magnetic field and residual eddy 
current, to be found. For example, the expression for the z- 
component of the magnetic field produced by a loop, combined 
from all four quadrants, is given by: 
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5 z (r) = J Bf(r) + 5;(r) + 5f / (r) ) 
2*f(r) = -^Xcos((2»-l)p)^ 

B\ (r) = ^ £ cos((2» - 1» Jdfc cos(fc)(AR, )3^_, (A) {© - p)/^, (kR, ) ( 1 0 > 

fi f ( r ) = ~^^Z(2«-l)cos((2n-l)^)sin((2/i~l^ }— cos(*Z(;t)) 

" =1 0 R p X 

{©(p - *)/ 2 , M (ip)^ (fa) + ®{x - /j)/^ (kx)K 2n _ x (kp)} 
where : 

= J^;(*)oo8[(2n-.l)/ p (z)]cos(fe), 

;f (id. 

3L-i(*) = J^;(^)cos[(2/2-l)/ 5 (z)]cos(fe) 

*f 

It will be recognized that all odd modes in cos((2h -\)<p) are 
present. These expressions hold for both isolated and 
communicating loops. For isolated coils 70, 80 that are closed 

10 on themselves, zf* ,z p / are the initial and final points of the 
primary (p) or shield (s) loop in the upper-half of the first 
quadrant. Communicating loops 72, 82 which connect with at 

least one coil jump 74 have zf* s as initial points on the 
primary (p) or shield (s) loop in the upper half of the first 
15 quadrant, and z p / s =L p ^ . From Equation (2) the result for the z- 
component of the residual eddy current is: 
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/,%,-. 00 = - fe, (*) sin(Az)cfe, 

gr i 1 

O)-- Or- W.^^^W-ttW^-^^W (12) 
sin[(2n — r p 

+ — 5 Hnf L ^*.*#.^«)«»(*»;)- : 5 L «(*.^.^«)«><fa;) 



+ (2«-l) 2 f— cos(AZ(x)) 1 2»-}( kx ) 1 v 



where the following notations are introduced: 



iO 2 



RJ I^ikRcsY 



(13) 



J^sin[(2 W -l)/ P(J (z)]sin(fe) 



The above approach can also be adapted to the design 
of a self -shielded axial gradient (Z-gradient) coil. In this 
case, it is sufficient to find the J 9 component of the current 
density on both the primary and the shield coils. This 
component does not necessarily vanish at the coil edges. In 
other words, the current continuity equation VJ = 0 does not 
imply that the J 9 component of the current density should be a 
continuous function of z even at the edges of the coils. 
Therefore, the values of the current densities at the edges of 
the coils can be treated as free parameters during the design. 
Thus , 
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J(r) = e^(r) 

J,<?) = {f/\z)S^R p )®{z+L p )®{L p -z) + f£\z)S(p-R 5 )®(z+L s )®(L s -*)} 

(14) 

where the functions f$ % *(z) are suitably represented as a Fourier 
5 series : 

ti P \z) = sin(^-z), f?\z) = 2>*> sin(^z) (15) 

with the coefficients a** 0 , a ( n s) being determined from 
10 extremization of the functional: 

^=*-£(*^ 

3 (16) 

(//(^)-^ p )a, -(/;(^)+4)A s 

where: E is the total energy of the gradient coil and the cold 
is shield; ^ co// (r.) is the z-component of the coil magnetic field at 
a set of points r. inside the DSV where the field has 
predetermined values B t ; .©^(r^is the z-component of the 
magnetic field produced by the cold shield eddy currents at a 
set of points r,on the surface of the DSV; B c z ° a (rj) is the z- 
20 component of the coil f s magnetic field at the same set of 
points; and C } is a set of numbers that govern the residual 

eddy current effect. Finally, JL t and Aj are the corresponding 

Lagrange multipliers. The quantities A p , A s are the desired 

values of the <|>-component of the current densities at the edges 
25 of the primary and shield coils, respectively. In Equation (16) 
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5 



10 



the shielding condition on the surface of the cold shield 20, 
namely: 

5 />( r )U^ = ° (17) 

has been incorporated. This leads to a continuous residual eddy 
current distribution given by: 



f? 00 = -)Fcs (*)sin(fe)<fe 

Fcs(k) = -Sh(k,R p ,R^F p ^ (18) 
F p Ak)= ]ff*{z)sm(kz)dz 9 Sh(k 9 R,R cs )= 



Optimization of the functional given in Equation (16) with 
respect to a[ p) and a< s) determines a suitable solution by 
adjusting the values of the current densities at the edge of 
the primary coil and the shield coil . After a solution is 

15 found the continuous current distributions on both primary and 
shield coils are discretized such that each loop on the primary 
coil carries the same current I, and each loop on the shield 
coil carries opposite current. The discretized coil is tested 
for its performance: gradient strength at isocenter, linearity 

20 and uniformity of the gradient strength over the DSV, slew 
rate, and residual eddy current effect. 

An exemplary transverse gradient coil design is 
described. The exemplary coil is designed for a magnet bore 
having a length of 1.2 meters and a bore radius of 0.45 meters. 

25 Design dimensions of the exemplary transverse gradient coil are 
as follows: R P = 0.34385m, L P = 0.434 meters, R s = 0.435224 
meters, L s = 0.555 meters. In all examples described herein the 
target eddy current effect was chosen to be -0.22%. 
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The coil is aligned such that the flare angle a shown 
in FIGURE 3 is 36.9°. In designing the X-gradient coil the 
field constraints shown in Table 1 were used- 



Table 1. 



p-position 
(m) 


CD- 
position 


z-position 
(m) 


B z (Tesla) 


0. 001 


0.000 


0.000 


2.7e-5 


0.180 


0.000 


0.000 


4.865285306501 
26e-3 


0.001 


0.000 


0.200 


2.069587125e-5 



The first constraint corresponds to the gradient strength of Gx 
= 27mT/m at the isocenter, the second constraint is assigned 
for no-rollover of the z-component of the gradient magnetic 
field, and the third constraint corresponds to -2 3% of non- 
uniformity of the gradient strength at z = 0.2 meters. The 
residual eddy current effect constraints were imposed at a set 
of points, which are given in Table 2. 



Table 2. 



z-position (m) 


p-position (m) 


0.0128 


0.2495 


0.06304 


0.2373 


0. 1296 


0.1904 


0. 19904 


0.0245 



The radius of the cold shield 20 was chosen to be R cs = 0.475/w . 
The coil jumps 74 between primary and shield coils 16 , 18 at 
the ends of the gradient assembly result in a significant 
reduction in stored energy and a change of residual eddy 
current effects. A study was done of total energy vs. the 
current density at the edge of the primary coil (quantity A in 
Equation 6) and the residual eddy current effect (quantity Cj , 

the percent measurement of the residual eddy current effect, in 
Equation 6) . 
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FIGURE 8A shows the plot for the energy of the 
system, rescaled to the gradient strength of 30mT/m and 
corrected back to account for the loss of the gradient strength 
due to the residual eddy current effect at the isocenter of the 
5 coil, as a function of the current density at the edge of the 
primary coil when the residual eddy current effect is equal to 
an exemplary -0.22%. FIGURE 8B, shows the plot for the energy 
of the system, rescaled to the gradient strength of 30mT/m and 
corrected back to account for the loss of the gradient strength 

10 due to the residual eddy current effect at the isocenter of the 
coil, as a function of the eddy current effect. In FIGURE 8B 
the current density at the edge of the primary coil is equal to 
an exemplary 17.0kA/m. FIGURE 8A shows that the total energy at 
the gradient strength of G x =3 0mT/m after accounting for the 

15 residual eddy current effect at the isocenter of the coil, has 

a minimum at a certain value of the quantity A for a given 

value of the residual eddy current effect. This value is 
suitably adjusted such that when the current densities on the 
primary and shield coils 16, 18 are discretized, the numbers of 

20 turns are integers. 

FIGURE 9 shows the results for the continuous current 
distributions shown in the unfolded coil layout of FIGURE 5. 
The continuous current distributions shown in FIGURE 9 were 
discretized with 29 turns on the primary coil 16 and 19 turns 

25 on the shield coil 18 per quadrant with 14 coil jumps 74 per 
half of each quadrant. Each closed loop carries current 
7 = 398.64^4 in order to produce gradient strength G x - 30.0mT/m 
at isocenter. The current layout per half quadrant is shown in 
FIGURE 5. Henceforth, P and S denote the whole primary coil 16 

3 0 and the whole shield coil 18, respectively, which together with 
the coil jumps 74 provide 50cmx40cm imaging volume. 

The discretized current was used to calculate 
magnetic field quality characteristics. FIGURE 10 shows 
X-gradient field non-unif ormity (taking the transverse gradient 

35 coil as producing an X-gradient; of course, the transverse 



WO 2004/046745 



PCT/IB2003/004802 



gradient coil can also be suitably oriented to produce a 
Y-gradient or other transverse gradient) , while FIGURE 11 shows 
X-gradient field non-linearity. 

FIGURE 12 shows the z-component of residual eddy 
5 current produced by the X-gradient coil having the current 
paths shown in FIGURE 5. The residual eddy current of FIGURE 12 
and the discretized loops shown in FIGURE 5 were used to 
calculate a residual eddy current effect over 50cmx40cm volume. 
The comparison of this result with that obtained from the 
10 continuous currents is shown in FIGURE 13. This figure shows 
that the discretization procedure does not introduce 
substantial error in terms of magnitude and its variation over 
the volume. 

With reference to FIGURES 14 and 15, the transverse 

is gradient coil whose layout is shown in FIGURE 5 is readily 
adapted to provide a discretely variable field of view. In a 
suitable modif ication, fourteen isolated primary coils 70 near 
the eye of the primary coil fingerprint are electrically 
decoupled from the remainder to define an interconnected coil 

20 shown in FIGURE 14 and designated as P2 herein. As shown in 
FIGURE 15, a remainder PI of the primary coil 16 remains 
electrically connected with the shield coil 18. For 
convenience, the original coupled primary coil 16 is designated 
as P herein. Hence, P=P1+P2. Similarly, the shield coil 18 is 

25 also designated by S herein. 

With reference to FIGURE 16, a second shield coil S2 
is designed to compensate for the fourteen electrically 
disconnected loops comprising the primary coil P2. A radial 
position of the second shield coil S2 is 0.415224m. The current 

30 layout on the second shield coil is shown in FIGURE 16. The 
current on the second shield S2 is 1=398. 64A. Combination of 
the coil shown in FIGURE 15 with the second shield coil shown 
in FIGURE 16, where the current flows in the same direction as 
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on the primary coil, provides a reduced DSV of 30cm. The non- 
uniformity and non-linearity of the gradient strength of two 
modes gradient coils are shown in FIGURES 17 and 18 , 
respectively. The main mode coil corresponding to the layout of 
5 FIGURE 5 provides 50cmX40cm imaging volume, while the coil of 
FIGURE 15 combined with the second shield S2 shown in FIGURE 16 
provides 30cm DSV. 

If a second power supply is provided, the gradient 
coil 16, 18 can be further modified to provide three discrete 
10 fields of view. The second shield coil S2 is replaced by an 
eight-turn shield coil S2 1 positioned at the same radial 
location. A base mode, which provides 50cmX40cm imaging volume, 
is generated by the gradient coil shown in FIGURE 5 driven at a 
current I . 

15 A second mode, which provides 3 0cm DSV, is generated 

by the gradient coil shown in FIGURE 15 driven with a current 
I, combined with the second shield 32 shown in FIGURE 19 driven 

with a current 1/2. Hence, two power supplies (or a current 
divider) are used. 
20 A third mode is generated by the gradient coil shown 

in FIGURE 15 operating concurrently with the inner fourteen 
loop coil P2 of FIGURE 14. The current in the primary coil P2 

is reversed relative to the current in the primary coil PI of 

FIGURE 15. The second shield S2 of FIGURE 19 is driven with a 

25 current -I. In this third mode, all coils are driven in series, 
and one power supply is used. This configuration provides 2 6cm 
DSV. 

When two coils are driven in series, a slew rate of 
the series-combined coils is given by: 

30 

SR = S* — - — (19) 
Li 
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10 



where: S is a sensitivity of the coil: S^y* with G and I 

being the gradient strength and the current f respectively; V is 
a output voltage of the power supply; R is a resistance of the 
coil including the cable resistance; and L is an inductance of 
the coil including the cable and filters inductance. 

With reference to FIGURE 20, when two coils are 
driven in parallel, the slew rate of the combined coil can be 
found by considering the coupling equations between two 
electric circuits, given by: 



L x ^ + RJ l+ M^ = V x , L 2 ^L + R 2 I 2+ M^ = V 2 (20) 



where L i9 R 19 V 1 and L 29 R 29 V 2 are the inductance, resistance, and 
output voltage of the first and the second circuits, 
15 respectively, and Mis the mutual inductance between these two 
coils. Then the slew rate of the combined coil is given by: 

SR -^{PA -S 2 MW X -R 1 I l ) + (S 2 L l -S X M){V 2 
A = L l L 2 -M 2 



20 where S l9 S 2 is the sensitivity of the first and second coils, 
respectively. 

Moreover, the gradient coils of FIGURE 5, including 
the primary coils 16 and secondary or shield coils 18, which 
are also designated herein as P and S, respectively, are 
25 readily modified to provide a continuous field of view. 

Starting with the initial coil system P+S, which 
provides the imaging volume with largest FoV, a variable FoV 
system can be achieved in the following way: (1) electrically 
disconnect several turns near the eye of the primary coil (this 
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defines the coil P2, for example as shown in FIGURE 14) and the 
remaining primary coil is PI = P - P2 (for example as shown in 
FIGURE 15) . The coil PI includes communicating coil turns 72 
that communicate with communicating coil turns 82 of the shield 
5 coil 18 by coil jumps 74 at the ends of the gradient. The 
positions of these interconnects are known and fixed. 
Optionally, the coil Pi also has some isolated turns that are 
only on the inner cylindrical surface 60. 

A second shield SI is designed to conform with the 

10 following conditions: (i) the number of conical jumps on this 
coil and their angular and axial positions match those on the 
primary coil PI; and (ii) shield SI together with the primary 
coil PI form a well -shielded pair. The system PI + SI provides 
a reduced DSV1 < DSV relative to the gradient coil 16, 18 of 

15 FIGURE 5. 

A second shield coil S2 ■ 1 is designed for the coil 
P2 . The second shield coil S2 ■ 1 together with the coil P2 form 
a well shielded pair. The combination of SI and S2 1 1 is 
equivalent to the initial coil S. The original coil P is not 
20 used (it was an intermediate process operation in the design 
procedure). Rather, two well-shielded pairs of coils: PI + SI 
and P2 + S2 1 1 are the final result of the design. 

These coil sets can be used as to provide either a 
discretely selectable FoV as in the previous examples or a 

25 continuously variable FoV. For a discretely selectable FoV, the 
coil (P1+P2)+(S1+S2 I ') provides the initial imaging volume DSV. 
The coil Pl+Sl provides DSV1 < DSV. The coil (Pl-P2)+ (S1-S2 11 ) 
provides DSV2 < DSV1 . 

With reference to FIGURE 21, for a continuously 

30 variable FoV, the coil sets (Pl+Sl) and (P2+S2 11 ) provide a 
continuously variable FoV when used in conjunction with two 
power supplies. In this arrangement, one coil set (Pl+Sl) is 
driven by a current I from a first power supply, while the coil 
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set (P2+S2 11 ) is driven by the current a*I from a second power 
supply. The parameter a adjusts the size of the FoV. 

With reference to FIGURE 22, another gradient coil 
with a continuously variable field of view is described. In 
5 this approach, a well-shielded coil is designed in accordance 
with process operations described previously, including current 
transfer between inner and outer cylindrical surfaces 
corresponding to primary and shield gradient coils, 
respectively, across a non-planar current-sharing surface. 

10 FIGURE 22 shows a suitable continuous current density for such 
a coil which produces a 32 cm DSV. 

A correction or shim gradient coil is designed in 
such a way that it produces zero gradient strength at an 
isocenter of the coil and together with the 32cm DSV gradient 

15 coil produces a 50cmx40cm DSV. By itself the correction coil is 
well shielded. 

The correction coil is suitably designed as follows: 
a well-shielded gradient coil is designed, including a 
non-planar current -sharing surface, that provides 50cmX40cm 

20 imaging volume and has the same gradient strength at the 
isocenter as that of the 32cm DSV coil. The value of the 
current density at the edge of the primary coil is the same as 
in the gradient coil that provides 32cm DSV. The difference 
(subtraction) between the current densities (on primaries and 

25 separately shield coils) of 50cmx40cm and 32cm DSV defines the 
correction coil. The values of the current density at the edge 
of the coils are adjusted in such a way that allows proper 
discretization. An example of the current density on the 
correction gradient coil is shown in FIGURE 23. 

30 The current distributions of the 32 cm DSV gradient 

coil are discretized to generate a current layout such as that 
shown in FIGURE 24. Similarly, the correction or shim coil 
current distributions are discretized to generate a current 
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layout such as the primary and shield shim coil current layouts 
shown in FIGURES 25 and 26, respectively. 

Variation of the amount of current in the correction gradient 
coil allows continuous variation of the FoV. The driving coil 
5 interconnections of FIGURE 21 can be used in driving the 32cm 
FOV gradient coils of FIGURE 24 and the correction coils of 
FIGURES 25 and 26. In this approach the coils PI and SI 
correspond to the well-shielded coil set of FIGURE 24 that 
provides the 32cm FoV, and P2 and S2' 1 correspond to the well- 

10 shielded correction coil pair shown in FIGURES 25 and 26. 

FIGURES 24-2 6 describe an arrangement in which the 
primary coil (FIGURE 24) produces a base field of view (32cm) 
that is continuously enlargeable by cooperative operation of 
the correction coils of FIGURES 25 and 26. Alternatively, a 

15 base coil can provide a large base field of view which is 
continuously reducible by cooperative operation of correction 
coils . 

Moreover, although design of transverse X-gradient 
coils are described herein, those skilled in the art recognize 

20 that the design of a transverse Y-gradient coil is 
substantively similar to design of an X-gradient coil, 
differing by a rotation of the whole gradient assembly by 90°. 

The gradient coil design procedures described herein 
are readily modified to accommodate other design aspects. For 

25 example, local forces on the gradient coil can be reduced by 
introducing an additional constraint that the net (thrust) 
force on the gradient coil be zero. This constraint 
incorporates information about the main magnet that the 
gradient coil will be used with. Reduction in the local forces 

30 on the gradient coil can reduce gradient tube deflection when 
the coil is operated. 

The invention has been described with reference' to 
the preferred embodiments. Obviously, modifications and 
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alterations will occur to others upon reading and understanding 
the preceding detailed description. It is intended that the 
invention be construed as including all such modifications and 
alterations insofar as they come within the scope of the 
5 appended claims or the equivalents thereof. 
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